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Subsidence of objects that rest on the seafloor is an 
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expected but not well understood phenomenon. Subsidence is 
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the settlement of fixed structures and mobile structures 
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(such as oceanographic packaces) that are placed on the 
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seafloor. A major concern involving bottom resting systems, 
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in relatively shallow water (less than 200 meters) is wave 
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induced motions. The system response of the structure in 
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this dynamic environment can create substantial loadings at 


the structure-seafloor interface. Past works have addressed 
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wave-induced and current scour as the major factors causing 
structural subsidence, neglecting the effects of pumping - 
erosion, 

The research described is directed towards indentifying 
and understanding the phenomenological aspects of pumping- 
erosion, The work describes a new hypothesis concerning 
liquefaction as being the precursor to subsidence in a two 


cycle pumping-erosion process. Due to the complexity of 
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structure-soil interaction a multivariable linear regression 


was performed on the data to substantiate experimental 
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CHAPTER | 
INTRODUCTION 


Structural subsidence can range from a few inches to 
several feet causing results varying from slight tilting and 
displacements of a structure to total and catastrophic 
structural failure. As a result of the petroleum and gas 
industries use of seafloor supported structures (i.e. jack- 
up rigs) much investigation has been directed at the 
structure-soil system and interaction. Although this area 
has been intensely investigated, the focus of attention has 


been on large oil rig structures. 


The use of relatively small mobile submerged structures 
for oceanographic and industrial purposes has increased in 
the past several years. The ability to successfully deploy, 
operate and retrieve these instruments and equipment 
packages is largely dependent upon the structure-seafloor 
interaction, Bottom resting systems are subject to the 
dynamic response of the structure in the ocean regime which 
can produce significant dynamic loadings at the seafloor. 
The stability of a submerged seafloor structure is dependent 
upon the subsidence and orientation of the structure while 
subject to wave and current forces. These forces produce a 


footing response in a variety of bearing, tension, sliding, 
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‘ rocking and torsional strains in the marine soil supporting 


the structure (Figure 1). 


The central issue of this report involves the 

< investigation of the phenomenon and mechanisms of pumping- 
erosion and its contribution to subsidence. Pumping-erosion 
is inherently non-linear as the structure is intermittently 
coupled and decoupled from the soil interface. In addition, 
the effects of the period of oscillation, amplitude of 
displacement, footing characteristics, hydrodynamic effects 
and soil parameters are simultaneously involved resulting in 


a very complicated and complex problem. 


The basic pumping-erosion cycle involves an uplifting 
force vertically displacing the footing off the seafloor. 
Subsequently, a downward force will dynamically load the 
structure onto the soil again. This over simplified cycle 
is a Trojan horse, as it involves several discrete and 


intertwined forces which are not readily apparent. 


Consideration must be given to the geometry and 
dynamics of the footing, The size, shape, and weight of the 
structure will determine such characteristics as center of 
mass, center of buoyancy, drag and lift forces and moments 
of inertia which will determine the structures response to 
the exciting force. The wave and current environment wil] 


produce the exciting force characteristics of period, 
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amplitude, and duration of excitation. Hydrodynamically, 
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effects of exit and entrance velocities, water density and 


¢ 
ata nats hae 
fa he ete 


ae, 


viscosity as well as the added mass of the interacting 


structure will influence the structure response and the 
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A sediment transport and settlement. 
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Soil parameters such as soil density, permeability, 
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r void ratio, grain size and shape, and pore water pressures 


will determine seafloor response to the structure-soil 
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interaction. 
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Traditionally subsidence was attributed to the effects 
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of scour resulting from the driving force of wave and 


1 


currents with little regard to pumping-erosion. Any nen 
inference to pumping-erosion usually surrounds the latter ae 
half cycle where the water is displaced or "pumped" from aay 
under the footing as it makes contact with the soil surface. “ie 
The soil transport has been attributed to the "jetting" nee 
action of the water escaping from under the footing. ed 
& 
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But what of the first half of the cycle of the footing 
when it is uplifted from the seafloor ? What forces are 
involved ? And what contributions do these forces make to 
the process of subsidence ? As the footing is uplifted, 
negative pore water pressures can be developed in the soil 


under the footing, possibly to the point of liquefaction. 
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Liquefaction, or fluidization, occurs when the pore- 


water pressure in a soil is equivalent to the total stress 
of the soil; in other words when the effective strength of 
the soil equals zero. Liquefaction of a soil requires 
considerable pore water pressures as compared to the 
relative ease of fluidization, which reduces the effective 
strength of the soil by the increase in the pore water 
pressure combined with the uplifting effect of the upward 
seepage flow. If the bed is considered deformable, the 
fluidization would then weaken the grain skeleton and 
distort the grain bed surface causing the grains to have a 
greater hydraulic profile making them readily available for 


transport during the downward half cycle of the footing. 


Thus, it is hypothesized that subsidence is a two phase 
cycle. The first phase occurs on the uplifting of the 
footing causing sufficient negative pore water pressures to 
fluidize the bed which weakens and distorts the bed making 
the upper grain layers more susceptible to sediment 
transport. The second phase, occurs on the downward stroke 
of the footing where the velocities developed under the 
footing transport the sediment. The transport to the 
sediment will be accentuated at the edges where the water is 
jetted out from under the footing. The following research 
is designed to investigate the hypothesis on a 


phenomenological basis. 
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The objectives of this research are as follows: 


1. To determine which geometric, 
hydrodynamic and kinematic parameters are involved 


in pumping-erosion subsidence, 


2. To evaluate the relative importance of 
the geometric, fluid and kinematic parameters 


involved in pumping-erosion subsidence. 


3. To determine if sufficient negative pore 
water pressures Can be developed under the footing 
to fluidize the foundation bed to the point of 


liquefaction. 
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CHAPTER II 
. LITERATURE REVIEW 
x 
e There are a limited number of references dealing with 
ia) 
soil/water/structure interaction. Even fewer address the 
‘ 
a soil/water/structure interactions and forces involved in eae 
7 = ~n 
- structural subsidence. Furthermore, there are only two hee 
* wets 
sited works that discuss the topic of pumping-erosion. Most er 
aa? 
Re structural subsidence is attributed to he effects of scour. om 
\ ad on 
Herbich et al.(3) has compiled an extensive biography on the tary 
< an 
he topic of conventional scour. nt 
aa 
yore *, 
oa ona 
Lad “ 
Reimnitz and Kempema (10) in 1981, describes the ean 
“ee 
° n f 
he effects of dynamic ice wallowing in the Alaskan nearshore *'e 


fi 
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areas. The report hypothesizes the occurence oftwo types of Nin 
a 
Kn hydraulic processes. In the first process the ice plays a “AN 
OL woo 
wih 
passive role, acting as a flow obstacle in the current and ey 
re wave regime, This particular aspect would account for the 
A, more conventional scour around structures. The second 
= 
hydraulic process involves ice playing an active role, 
Me either by simple vertical oscillations or by wallowing 
wh 
(rocking). Dahlberg (1) discusses pumping erosion as it 
i pertains to structures. He states that pumping erosion is 
wt 
i associated with excess pore-water pressures set-up in the wns 
Pa 
. foundation soil during storm periods. Dalhberg suggests ae 
oe 
Aa) ras 
" “* 
o 
~ 


De rae? te Le ae ale ‘ee CP iid a Panta hs! af att ne”, Pied al »,! te Ag ye Mra SO ao 


eae ie. 


po 


. ‘w a. 


that excess pore-water pressures develop to balance the 
overturning moments on the base of the structure. It is 
further stated that if there is free communication between 
the foundation soil and the seabed, the pore-water gradient 


may be high enough to liquefy the soil locally. 


It was verified that a North Sea platform, Frigg CDP1, 
which had no protective scour skirts, was subject to 
pumping-erosion. Divers observed periodic puffs of sediment 
around the periphery of the structure base when the 


sediments were carried in suspension as the structure rocked 


in the storm environment. It is noted that the pumping- 


erosion described was solely attributed to the positive 
pore-water pressures escaping above the mudline due to the 


pressure gradient. 


An extreme case of pumping-erosion involved the 
Christchurch Bay Tower located in 8.4 meters of water. The 
platform was subject to severe storm conditions. Due to the 
absence of protective scour skirts the platform was 
undermined by scour, which lead to a free rocking motion of 
the structure, resulting in pumping-erosion on the 


foundation soil. 


Pumping-erosion of gravity structures is now rare,due 


to their size and the use of scour skirts. Most settlement 
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Teramoto et al. (13) discussed the scour encountered w..h es 
various sit-on-bottom type of structures. The study 
investigates the various scour patterns and contours for i 


various footing designs and configurations. x 
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However, pumping erosion is more applicable to 


; lightweight mobile structures (i.e. oceanographic research ye 
: instrumentation packages which can be deployed from research 2 
5 vessels) and is the focus of this topic. , 
Several works have dealt with wave induced pore water Z 
- pressures, Approaches have varied from considering the soil ‘ 
skeleton to be rigid, to being compressible and pore water 
as being compressible or incompressible. Governing e 
s equations were usually derived from Darcy's Law of flow : 
i through a permeable bed or Biot's Law for three dimensional * 
consolidation for a poro-elastic material. In addition, > 
; various studies have evaluated whether residual pore-water ; 
; pressures exist in cohesionless soils or whether the pore oy 


pressures are merely transient. 


Oldinziel and Brink (9), indicate that an upward flow oe 
of water through a porous sand bed reduces the apparent a 


weight of the sand particles and therefore reduces the sand 
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particles’ stability. The results of the study conclude oS 
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that upward flow (blowing) through the bed increases the 


rate of sediment transport. The study relates the upward read 
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Martin and Aral(5), 1971, indicates that the seepage 


force on a surface grain is only 50% of that for an embedded bd 
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whereas the downward seepage will increase stability. 
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CHAPTER III 
PROCEDURES AND INSTRUMENTATION 


When non-fixed seafloor systems are subjected to wave 
and current forces of sufficient levels, the structure will 
oscillate. The wave pressure forces, in conjunction with 
the current forces will lift the structure and then will 
load it again on the seafloor. The response of the 
structure and subsequen. subsidence are a function of 


several geometric, hydrodynamic and kinematic factors. 


Due to the amount of limited data concerning pumping- 
erosion subsidence a simple plate type footing employing 
purely vertical oscillations was chosen. The tests were 
conducted in a steel tank with a glass observation window. 
The dimensions of the tank were 8 feet long, 4 feet wide and 
4 feet deep. The glass observation window was 57 inches 
long and 3/4 inches thick. A sand bed, 11.5 inches deep was 
placed in the steel tank. The test section (4 feet by 4 
feet) was 11.5 inches deep. The remainder of the tank had 


9.5 inches of sand overlying a false bottom (Figure 2). 


The sand utilized for the test was an Ottawa sand with 
é specific gravity of 2.65 and a mean grain diameter (d50) 


of .22mm (0.087 inches) (Figure 3). The sand had a static 


(pluivated) submerged angle of repose of 16.50. 
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Two steel plates were utilized in the experiments. A ror 
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2 inch square plate (1 square foot), 3/4 inch, 25 pound 
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plate was utilized for experiments 1 through 6. A 10 gauge aon 
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at 
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(0.132 inch) thick, 16.97 inch square (2 square feet) plate pace 
attached to the 1 square foot plate was utilized for the 


experiments 7 through 12, 


The plate was uplifted by a cable and pulley system near 
Sec 


coupled to a 1/2 horsepower D.C. variable speed motor via a : a 
rotating cam follower that engaged lever arm attached to the i 
cable (Figure 4). This particular arrangement allowed the Aes 
plate to be uplifted from the sandbed to a predetermined ty 
amplitude and then the plate was released and allowed to ea 


free fall onto the test bed. The next cycle would engage 
the cam follower on the lever arm lifting the plate again, 
repeating the cycle. The plate was thus able to follow the aN 


contours of the scour caused by the pumping-erosion, 


The amplitude of oscillation (1 or 0.8 inches) was 
obtained by eccentrically locating the cam follower from the 
center of rotation of the motor shaft. The selection of the 


amplitudes was restricted by equipment limitations. 


Oscillation periods cf 5, 10 and 20 seconds were 
selected to simulate general storm wave conditions that 


might be experienced by seafloor systems, 
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Pore-water pressures in the test bed were measured by a 


Statham P-22S-350 bi-directional pressure transducer (2 


psid). The transducer was placed outside of the tank at the pee 
oN 4 


same equivalent hydrostatic depth as the sand bed. The ae 
tet 

transducer detected pressure fluctuations through the active “ES 
aa w5 

port via a 1/8 inch clear PVC tubing coupled to a needle may 
Poe 


shaped copper probe fitted with a # 200 mesh screen to 


ay 
restrict the entrance of sand particles into the transducer oe 

-. 
port. The transducer probe was centered under the plate- re 
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footing at a depth of 0.276 inches (7mm). 


To insure proper measurements, the active transducer 
port, PVC tubing, brass coupling and sensing probe were 
thoroughly cleaned with a strong detergent to remove any 
dirt, grease or oil that may have accumulated. After 
rinsing, the active pressure port and extensions were 
carefully flushed and filled with de-aired water. This 
procedure was followed to eliminate the possibility of 


entrapped air within the pressure sensing system, 


A static load test was utilized to calibrate the 
pressure transducer. Readings were taken with the 
transducer submerged in still water at a predetermined 
depth. At a second predetermined depth readings were taken 
again, Since the pressure relationship is a linear, an 


average calibration factor of %.043 psi/mv was obtained. 
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The 5 volt D.C. excitation voltage was supplied by a 
HP17403A pre-carrier amplifier. The output was recorded by 


a HP7402A dual channel strip chart recorder. 


To record the near bed shear velocities at the plate 
edge, a Thermo-Systems Inc (TSI), hot film anemometer system 
was employed. The voltage output of the anemometer is 
proportional to heat loss across the cylindrical filament 
(film) as a result of the fluid flow across the filament. 
The actual velocities were determined from calibration 


curves developed prior to each experiment. 


The calibration was performed by attaching the hot film 
probe to a motorized trolley carriage which rode on a rail 
system above a 120 foot wave flume. The carriage was 
operated at several different speeds over a predetermined 
distance. The voltage output from the TSI anemometer 
monitor was displayed on a Beckman 200 digital voltmeter and 
recorded, A calibration curve of millivolts versus ft/sec 
was thus developed. Temperature corrections based on 
Reynolds Number similarities were employed to compensate for 
the difference in water temperatures between the wave flume 


and the test tank. 
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The upward plate velocities (VPU) and the downward 
plate velocities (VPD) were determined from the hot-film 


anemometer velocity record. 


The scour depth readings were recorded manually by 
observing a marker attached to a stanchion affixed to the 
plate traveling across a stationary rule. Horizontal 
deflections were minimized by extending a guide tube over 


the square tubing shaft welded to the plate surface. 


Twelve experiments were conducted at various 
combinations of period (T), amplitude of oscillation (A) 


? 


and plate dimensions (area (AP); weight (W); and thickness 


(T) ). 
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CHAPTER IV 
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EXPERIMENTAL RESULTS 
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A total of 14 experiments were conducted. Two 
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preliminary experiments were performed to observe the ve 
effects of the transducer probe position under the plate. se 


As one might expect, the highest negative pore-water 
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pressures were recorded at the center of the plate. 
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Therefore, the center of the plate was chosen for the probe 
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placement. It was the goal of these experiments to obtain 
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a general understanding of the mechanisms involved in the 
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phenomenon of pumping-erosion. The experiments were also 
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designed to investigate whether negative pore-water ~ 
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pressures and possible fluidization of the foundation bed 
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played a significant role in pumping-erosion. Therefore, it 
was necessary to discern the effects of the experimental 
parameters in pumping-erosion subsidence. Following is a 


table listing the controlled variables for each test run: 
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TABLE 1 - Controlled Parameters i 
sce eects wed a ale a te en : . 
~\ 
EXP PERIOD AMPLITUDE PLATE PLATE PLATE " 
NO. (T) (A) AREA WEIGHT THICKNESS bs 
(AP) (W) (TH) 
sec - inches ft2 Tbs inches 
l 10 1 1 25 15 
2 20 1 1 25 75 eS 
3 5 1 1 25 75 
4 20 l l 25 75 ONG 
one 
5 10 1 1 25 75 twa 
as 
6 5 1 l 25 75 es 
ee 
7 20 2 80 36 132 SR 
8 10 2 80 36 Sa 
132 stn 
9 5 2 80 36 132 
10. 20 2 80 36 132 2 
it... 10 2 - 80 36 132 en 
i, 2 - 80 36 132 acs 
erate 
BN ei eee te ae Eee ee eS Hare en Te By 7, gt Sat 
xs 
The soil used was an Ottawa sand with a mean grain met 
diameter of .22mm, The average bed density of the note 
foundation bed (the area directly under the plate) was 3.30 a SS 
slugs/ft3. In addition the controlled parameters listed in eege 
pea 
Table 1, the following variables were also recorded: number see 
Ce tie he 
wt ts 
of cycles(CYCS), positive pore water pressures (PPWP), Pde 
; on 
negative pore water pressures (PPWS), plate exit velocities Sordi de 
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Figures 1 through 12 in Appendix A plot the subsidence te 
~e 
| versus time and cycles. Figures 13 through 24 plot the aioe 
subsidence versus the log of time and cycles. aon 
Figures 25, 26 and 27 are group plots of subsidence ne 
SOREN 
versus cycles and Figures 28, 29 and 30 are group plots of wes 
ao - 
subsidence versus time for the periods of 5, 10 and 20 os 
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seconds respectively. Figure 31, plots the subsidence 
curves versus cycles (Figure 32 versus time) for the first 
six experiments (1 sf plate). Figure 33, plots the ie 


subsidence versus cycles (Figure 34 versus time) for the 
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x last six experiments (2 sf plate). 


Figures 5A through 5L show typical strip chart output - 


records for the 12 experiments. The pore water pressure 
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- record output is recorded on channel 1 and the hot film ee 
ef anemometer output record on channel 2. jee 
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CHAPTER V 
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DISCUSSION OF RESULTS 


As one reviews the results depicted by the figures 
(Appendix A), certain trends appear. The first observation 
is that the experiments that produce the most rapid 
subsidence have the larger plate areas, greater weights, 


larger amplitudes of oscillation, and shorter periods. 


. 


<W 


Secondly, the experiments that subside more rapidly also had 


higher negative pore water pressures and thus higher 
liquefaction indexes as derived by Sleath (12), where the 


point of fluidization is regarded as : 


water density 


ted 


bed density 


Pa a 


pore-water pressure through the transducer depth dy 


ve 
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Failure of the bed could also result from too large a 


horizontal gradient : 


. = 
y Bs 
be: Spy = peer: tani 8 ae 
“rg cy 0 a 
: where 3 = internal angle of friction 
“AS the plate is uplifted, as depicted in Figure 6, ; 6 
negative pore water pressures are developed locally in the pee 
Fk foundation soil directly under the plate. The data indicate Red 
. that the bed is indeed fluidized. Assuming that the pore- -. 
‘ water pressure is incompressible and that the grain skeleton ay 
m is compressible and deformable, the action of local aM 
fluidization of the bed would than cause distortion of the as 
. upper grain layers (see Figure 7). a 
~ AS 
2o8 
Some of the loose grains were observed to be carried oe 
e into suspension by the turbulent flow of fluid rushing in ny 
: under the plate in response to the local negative pressures. aa 
” The strip chart record reflects the velocity response to the ek 
: pressure fluctuations which are in agreement with the 
i Bernoulli equation 
: pz-}oV2- gx + 
a 


COPA EN Can a ae ee SE oa OF ‘eet : eae PO ACERT te a Pera nN En 
fe ‘a ‘e e a 5 re A 5 ie ts 7 i *' f. 5 


yy 
pata’ é fe 


Sat dat) tA Bel tat ctel Sy, 


bi emeth 


ov 
Pédd’. 


nearer ae Seae S 


Fre RT) 


VAS SaS 


Ca are 


“« 
“a 


ORM .- 
TION) JLT Cavmdy) 


BUNSS3dd VAL eHaN 


fe a 
Cie CAP POCO L 
Dass va @ 


a? 
“oe 


4, eo," 


en Oe ae | eo ; nee 
ae é ade Os one, PS Fis ete 4 aS ‘e a aad fe . ? . 


/ RNs os 
NOLITANS NIVY9 CaLYOLSI(] bee 


37 


Tr we wey 


‘ La 


( : 
2 
sok 
CX 4 
Cae: 

0 
eS 
Ox 


wr 


Site 
? 


— é 

ee 

SUIAVT ‘ 
NIVYD ee a 


GjINOLSIA MO14 NOILWZIGINT ‘" 


Se ee eo 


Bae tate ne tele 
. 
ate 


A% 


OE rie “ eps aD . Bee Fhe fe, 46, ; fee err Bek bt on oS o 


RRA 


o * : = . » : v e 
Se ee el te eee i Be SN NO} BR SS we ae ag MISES eee Ime aes ee 


During the downward cycle of the plate (Figure 8), the 
loosened grains as well as the sediment in suspension (as a 
result of the uplifting of the plate) can then be easily 
transported (Figure 8) at a much lower shear velocities (u*) 
then would be required for normal initiation as predicted by 
the Shields' diagram (6). An interesting observation, 
derived from the pressure and velocity records, is that the 
slower periods have more pressure fluctuations and thus more 
velocity fluctuations but with lower negative pore water 
pressures than do the shorter period cycles which have 
greater negative pore water pressures but with a faster 


change in pressure and velocity with fewer fluctuations. 


By careful observation three interesting details were 
discovered. First, a small ridge of sand developed along 
and under the perimeter edge of the plate. The second 
observation is that a trough was developed immediately 
adjacent to the perimeter edge of the plate (Figure 9). 
Thirdly, slope instability and sloughing was observed during 
the upstroke of the plate. These seemingly minor details 


may play a significant role in the pumping erosion process. 


The small ridge under the perimeter of the plate 
probably results from the outflow of water and transported 
material attempting to exit from the underside of the alate 
during the last fraction of a second before plate impact. 


This small ridge of sand may account for the large 1eqative 
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pore water spikes at the end of each cycle (Figure 10) in 
one of two ways. First, if the foundation bed is considered 
to be elastic, the negative pore water pressure may be due 
to plate-soil bounce with the small ridge acting to seal the 
plate so that the pressure response is localized. The 
second possibility, is that as the plate impacts on the 
foundation soil, the plate undergoes slight flexure and 
causes the negative pressure excursion. Again, the sand 
ridge would tend to seal off the plate to cause the event to 


be a local effect. 


The trough is caused by the jetting of water from under 
the plate (Figure 9). During the uplifting of the plate, 
the lower portion of the slope became unstable and sloughed 
onto the aforementioned trough. The entire slope then 
cascaded down to form a new equilibrium slope. The slope 
failure was probably due to the pressure/velocity field near 
the edge of the plate. As the water rushes in, responding 
to the negative pressures under the plate, sand grains at 
the base of the slope could be enveloped by the turbulent 
flow causing a slope failure. An alternate explanation, is 
that there is flow out of the slope in response to the 
pressure fluctuation, which would also cause the slope 


failure. 
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: In any event, the sand from the slope failure, which sete 
” should be large grains due to preferential transportation we 
r (winnowing), loosely fall into the trough only to be =a 
expelled and transported during the downward stroke of the wee 
plate. ee 
Topographical plots of the sand contours for the 12 EN 
- experiments are in Appendix C. Of particular interest is nan 
. the dynamic angle of repose of the scour mound slopes as aN 
listed in Table 2. aS 
: TABLE 2 DYNAMIC ANGLES OF REPOSE ee 
] 21 
‘i 2 24.3 
j 3 24.9 
4 24.2 
: 9 23.4 
a 6 25.8 
; / 2353 
a 
(i 8 22.2 
9 23.8 
10 2252 : 
; ll 22.9 ee 


me Le 222) 


Figure ll, ranks the experiments by the average initial 
liquefaction index. The ranking is in accordance with the 
observed subsidence rates (i.e. the tests with more rapid 
subsidence have higher liquefaction indexes). A 
relationship between the period, amplitude, plate weight and 
plate area ( A2Ap/WT2 ) was then developed and ranked 
(Figure 12). Curiously enough, the ranking was the same as 
the liquefaction index except for the positions of 


experiments 2 and 4, which were reversed. 


Although some general trends were developed from the 
plotted data and figures, a more sophisticated numerical 
analysis is required to ascertain the relationships between 


the experimental parameters. 


MULTIVARIATE PROCEDURES 


To determine the relationship and relative importance 
of the variables the Statistical Analysis System, or SAS, 
was employed. The SAS techniques were used to analyze the 
data because of the statistical integrity and ease of use. 
All of the statistical methods employed are explained in the 
SAS User's Manual : Statistics (ll). The data from all 12 
experiments were compiled into 234 observations, using 15 
variables. The dependent variable was subsidence (SUBSIO). 
The 14 independent variables were : neriod (T), amplitude of 
oscillation (A), plate area (AP), plate weight (W), plate 


“thickness (TH), bed density (BRHO), kinematic viscosity 
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(VISC), number of cycles (CYCS), negative pore water ay 


pressures (PPWS), liquefaction index (LI), positive pore- 


water pressures (PPWP), plate upward velocity (VPU), plate 2 
* 
downward velocity (VPD) and plate exit velocity (VE). Table Tere 
Ree 


3 summarizes the range of values for each of the variables. 
A list of all the variable values can be found in Appendix 


B. 


The SAS procedure, PROC CORR, was first used to 
determine simple linear correlations between the variables 
in the entire data record and are listed in Table 4. Due to 
the high correlation between the area of the plate (AP) and 


the plate weight (W) and plate thickness (TH), W and TH were = 


deleted form further analysis. Also LI was deleted from i 
ee key 

analysis due to the high cross correlation with PPwS. TOSS 
ate 

Furthermore, PPWP was not considered for further analysis ree 
ra ie 
eta 


because of its low correlation (-0.016) with subsidence. 


- 8 
a 


This is an interesting event as data suggests that positive 


we wt 


pore-water pressures have little to do with subsidence. 


The next step was to determine the strength of the 
Tinear relationship determined by computing R2 , the 
square of the multiple correlation coefficient (R-SQUARE in 
SAS). R2 also known as the coefficient of determination, is 
interpreted as the portion of the variability that has been 
accounted for by the regression formula. The closer R2 was 


to 1, the better the model equation fit the data. However, 
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a large R2 did not necessarily mean that the proper model 
was chosen for the data. Therefore, R2 should be used with 
caution, since it is always possible to make R2 large by the 
addition of more variables to the model. The R-SQUARE 
routine was configured to give R2 and C(p) values for all 
possible combinations of 3, 4 and 5 variables. The results 
are listed in Appendix D. C(p), Mallow's statistic, is a 
criteron related to the mean square of a fitted value (7). 
As C(p) attempts to indicates the bias of the prediction, a 
small] C(p) value is generally sought. The “best fit" model 
is determined by when C(p) first approaches p on a plot of 


C{p) versus p. 


Next, the remaining 19 variables were incorporated into 
Av ore ytine called STEPWISE. STEPWISE is a model building 
vam ars & oneonses variables that are to be included into 

Ttea'ly, a variable selection process should 
r any jtven nymber of independent variables, a 

av ove stuatica of predicted values produce the 
mm o* she squares, In PROC STEPWISE 
o4t ws gre aydted or qeleted one at a time, depending on 
tear t or 4a ewart .otion selection, until all 
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The STEPWISE procedure, using the backward elimination 


method, gave results that are suspect. 


It indicated that 


the best model was produced when all the variables were 


significant at the P < 0.001 level when the variable AP 


(plate area) was removed. 


This result is contradictory to 


the forward selection process which adds variables in the 


order of there significance. 


ranked the variables in order of significance in the 


following manner: 


TABLE 5 


Variables 
cYcs 

VPD 

VISC 

AP 


C 
185.049 
28.146 
20.663 
16.634 
14.272 
11.090 
9.307 


FORWARD SELECTION RANKING 


F 
381.754 
143.304 
8.842 
5.737 
4,209 
o<091 
3.2761 
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The significance of the ranking is apparent in the 
manner in which the forward selection was constructed, in 
that it added the variables in the following order : CYCS, 
VPD, VISC, AP, T, PPWS, A and VPU. Each additional) variable 


improved the overall fit of the model. 


The strength of a linear relationship, an indication of 
how well the model fits the data, is indicated by the R2 
value. The F-value is used by SAS to test the null 
hypothesis, which tests the linear association between Y and 
X1, X2, ... Xk . If the F-value is not significant (PROB F 
= 0.0001), for a system an "n" degrees of freedom (OF), the 
null hypothesis is accepted indicating that the model does 
not explain a significant portion of the variability. That 
is to say, as the F-value increases a greater portion of the 


variability is explained by the model (i.e. the spread of 


the model is lessened). 


The PROB >F indicates the significance of the model 
(or the strength of the variables as in the R2 output). The 
Tower the value of PROB>F , the more significant the model 


A p-value of 0.0001 is highly significant. 


The strength of the individual variables in indicated 


by, T and PROB>T, which denote the same inference as F and 


PROB>F, 
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Based on the above results, PROC REG was performed on 
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the 8 remaining independent variables. Variables A 


z, 
1 
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rig 


(amplitude) and VPU (plate upward velocity) were then 


vara eat 


removed and PROC REG was run again. The results are shown 


ar 


in Tables 6 and 7. 


Two variables, Fl and F2 were then created from the 


combination of the 8 variables where : 
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F2 = {(1/2VPD2 A W )/(144 VISC2 g2 T2 )}1/2 
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F3 and F4 were created from the log of Fl and F2 taking 
into account the number of cycles. PROC CORR, PROC RSQUARE 
and PROC REG were then utilized on Fl, F2 F3, and F4. The 


results are listed in Tables 8, 9 and 10. 


Plots of the models (Figures 13, 14 & 15) listed in 


Tables 6, 7 and 10, graph the predicted values from the 


models against the actual values. 
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CHAPTER VI 


SUMMARY AND CONCLUSIONS 


The main objective of this research was to determine 
which mechanisms were involved in the phenomenon of pumping- 
erosion, The secondary objectives, involved the evaluation 
of the relative importance of each of the 14 independent 
variables, and in particular, the significance of negative 


pore-water pressures. 


Of the original 14 independent variables, 6 variables 
(cycles (CYCS), downward plate velocity (VPD), viscosity 


(VISC), plate area (AP), negative pore water pressure 


LP ORD 


(PPWS), and period (T)) were considered significant in the 


wv 
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model constructed utilizing SAS. The number of cycles 
(CYCS) and the downward velocity of the plate (VPD) were 
highly significant and is in agreement with the accepted 


understanding of pumping-erosion and scour. 


The result of negative pore water pressure (PPWS) 
remaining significant throughout the analysis supports the 
hypothesis that negative pore water pressures are a 
contributing factor in pumping-erosion. This is 


understandable if the experimental results of the 


Liquefaction Index (fluidization index) associated with the 


negative pore water pressures is accepted. 


This leads to the conclusion that the fluidization of 
the foundation soil, due to the negative pore water 
pressures developed during the upward cycle of the plate- 
footing is a precursor to the scour developed during the 
downward cycle of of the plate-footing and significantly 


contributes to the overall] pumping-erosion subsidence . 


As with many initial investigations, this research 


poses more questions than it answers. Such as 


1) If negative pore-water pressures and the 
subsequent fluidization are significant to the 
phenomenon of pumping-erosion, then what is the 
distribution of the negative pore water 


pressures under the footing ? 


2) What is the distribution of velocities under 


the plate footing ? 


3) Can the footing be designed to alter the 


pressure and flow nets under the footing to 


reduce subsidence ? 
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4) What effect does the shape of the scour mound 


have on the negative pore water pressures ? 


5) How would the results vary in other 


cohesiontess and cohesive soils ? 


6) What are the effects on pivotal oscillations 


as opposed to vertical oscillation ? 


RECOMMENDATIONS 


To obtain a better understanding of pumping-erosicn, 
more research is needed in this area. If the results are 
valid, the effect of negative pore water pressures may 
affect all dynamic structure-marine soil interactions 


including submarine pipelines. 
For purposes of further research the following 


recommendations are proposed 


1) Employ hydraulic or pneumatic rams to load the 
soi] at a constant rate with a constant 


amplitude. 
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formation of the scour mound slope. 


Conduct the experiments with different grain 
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sizes of cohesionless soils. 


Conduct the experiments with cohesive soils. 
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Alvin Eugene Grimmig Jr. was born on May 24, 1951 in 
“Miami Beach, Florida of Alvin and Nora Grimmig Sr. He 
married Miss Debra Kay Rademacher on May 6, 1972. He 
enlisted in the Navy in June 1972 and received a commission 
in Civil Engineer Corps in 1978, He received a 
baccalaureate degree in ocean engineering from the 
University of Washington located in Seattle, Washington in 


December of 1978. 


Permanent address: 


LT Alvin E. Grimmig Jr. 
c/o Barbra Rademacher 

555 Filmore Ave. Apt.101 
Cape Canaveral, FL 32950 
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